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In this paper, based on the nonbinary graph state, we present a systematic way of constructing 
good non-binary quantum codes, both additive and nonadditive, for systems with integer dimen- 
sions. With the help of computer search, which results in many interesting codes including some 
nonadditive codes meeting the Singleton bounds, we are able to construct explicitly four families of 
optimal codes, namely, [[6, 2,3]]p, [[7, 3, 3]]p, [[8, 2,4]]p and [[8,4, 3]]p for any odd dimension p and 
a family of nonadditive code ((5,p, 3))p for arbitrary p > 3. In the case of composite numbers as 
dimensions, we also construct a family of stabilizer codes ((6,2 •p^,3))2p for odd p, whose coding 
subspace is not of a dimension that is a power of the dimension of the physical subsystem. 



I. INTRODUCTION 

Noises are inevitable and they cause errors in quantum 
informational processes. One active way of dealing with 
errors is provided by the quantum error-correcting codes 
(QECCs) [H El [3 H], which have found many apphca- 
tions in quantum computations and quantum communi- 
cations, such as the fault-tolerant quantum computation 
[5], the quantum key distributions and the entangle- 
ment purification ^|8j. Roughly speaking, a QECC is a 
subspace of the Hilbert space of a system of many physi- 
cal subsystems with the property that the quantum data 
encoded in this subspace can be recovered faithfully, even 
though a certain number of physical subsystems may suf- 
fer arbitrary errors, by suitable syndrome measurements 
followed by corresponding unitary transformations. 

An important family of QECCs is the stabilizer code 
[Hldnilll], which is specified by the joint -1-1 eigenspace 
of a stabilizer, an Abelian group of tensor products of 
Pauli operators. The stabilizer formalism has also been 
established in the nonbinary case [12l [13l [HI |T5j [16] and 
many good codes for systems of a prime or a power of 
prime dimension have been constructed [15| 1161 1171 [TBj , 
including the well-known perfect code with five registers 
[T7] for all dimensions. Though the majority of QECCs 
constructed so far are stabilizer codes, including the CSS 
codes [H], the topological codes [3D], color codes [3T], 
and also the recently introduced entanglement-assisted 
codes [22], there are a few exceptions called as nonaddi- 
tive codes ^M^M- 

The nonadditive code does not admit a stabilizer struc- 
ture and it should not be a subcode of some larger sta- 
bilizer code with the same distance otherwise it will 
be a trivial nonadditive code. Though difficult to con- 
struct and identify the nonadditive codes promise a larger 
coding subspace since less structured than the stabi- 
lizer codes. For qubits the nonadditive error-correcting 
code that outperforms the stabilizer codes has been con- 
structed [26] based on the binary graph states. A graph- 
ical approach [57], as well as a codeword stabilized code 



approach [28j , to the construction of binary additive and 
nonadditive codes has been developed based on the bi- 
nary graph states. 



The graph states [Ml [30] are useful multipartite entan- 
gled states that are essential resources for the one-way 
computing [31] and can be experimentally demonstrated 
[32] . The binary graph state proves to an extremely effec- 
tive tool [m [301 [33] in the construction of QECCs. The 
nonbinary graph states were introduced first in [30J and 
discussed in details in the case of systems of an odd prime 
dimension [3S1. It is also investigated in the context of 
universal quantum computation [36] and other applica- 
tions [37]. Recently an approach to construct QECCs 
based on nonbinary graph states has been introduced in 
[35] and some new codes are found via computer search 
for qubits and qutrits. 



Here we shall generalize the graphical construction of 
QECCs [27| to the nonbinary case based on the nonbi- 
nary graph states from which some analytic constructions 
are attainable. In Sec. II the nonbinary graph states are 
introduced. In Sec. Ill we introduce the concept of cod- 
ing clique for a weighted graph and show how it is re- 
lated to the construction of both stabilizer and nonad- 
ditive QECCs. In Sec. IV we present some codes found 
via numerical searches and the graphical versions of some 
known codes as illustrations. In Sec.V we construct ana- 
lytically four families of optimal stabilizer codes that sat- 
urate the quantum Singleton bound for any odd dimen- 
sion as well as a family of nonadditive codes ((5,p, 3))p 
for all p > 3. In Sec. VI we investigate the graphical codes 
arising from composite systems and construct a family of 
stabilizer codes ((6,2 •p^,3))2p with p being odd. The 
graph states of systems composed of coprimed subsys- 
tems are in a one-to-one correspondence with the direct 
product of the graph states of subsystems. 
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II. NONBINARY GRAPH STATES 

Here we shall consider a general system with p levels, 
a qupit for short, where p is arbitrary. We denote by 
Zp = {0, 1 . . . ,p — 1} the ring with addition modulo p. 
Under the computational basis {\i)\i G Zj,} of a qupit, 
the generalized bit shift and phase shift operators read 



(1) 

Obviously XP = ZP = J and ZX = luXZ. Here we 
denote by X the Hermitian conjugate of X. The compu- 
tational basis |/) is the eigenstate of Z with eigenvalue 
while the eigenstate of X with eigenvalue reads 



1^. 
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(2) 



A Zp-weighted graph G = {V, T) is composed of a set 
V of n vertices and a set of weighted edges specified by 
the adjacency matrix T € Zp^", an n x n matrix with 
zero diagonal entries and the matrix element Tab G Zp 
denoting the weight of the edge connecting vertices a 
and b. The graph state associated with a given weighted 
graph G = {V,T) of a system of n qupits labeled with V 
reads EHl 



(3) 



a,beV 



Here we have denoted by Z^ the set of all the vectors 
s ~ (si, S2, . . . , Sn) with rt components Sa G Zp (a € y), 
by |s) the common eigentsate of all the phase shifts Za 
(a £ V) with eigenvalue , by 



(4) 



the joint +1 eigenstate of all bit shifts Xa {a E V), and 

by 



Uab= ^'^N)(*|a® |i)(i|fc 



(5) 



the non-binary controlled phase gate between two qupits 
a and b. The non-binary graph state |r) is also the unique 
(up to a global phase factor) joint +1 eigenstate of the 
following n vertex stabilizers 



ga = XaY[{Zbf-\ aeV. 



(6) 



For s e Z^ we denote X" = X'^X^^ ■ ■ ■ and simi- 
larly for the phase shift operator Z®. Obviously 



gs= UiGaY- = uj'-''-^-'x'z' 



(7) 




FIG. 1: Graph, graph states, and graph-state basis, (a-c) 
Star graph on 3 qubits and graph-state basis; (d) Star graph 
on 4 vertices; (d-e) Loop graph on 5 vertices. Colored vertices 
represent also a vector in Zpi white, black, blue, and red 
vertices having weights 0, 1, 2,p — 1 respectively. Black edges 
have weight 1. 



is also a stabilizer of the graph state for arbitrary s g Z J , 
i.e., Gs\T) = |F). All the stabilizers of the graph state 
belong to the generalized Pauli group for qupits 

p„ = {e-*5(p-i)"-*A'^Z*|s,t e Z^} X e Zp}. (8) 

The graph-state basis of the n-qupit Hilbert space 7i„ 
refers to {|Fc) = Z'=|F)|c e Z^}. Under the computa- 
tional basis a graph-state basis looks like 



(9) 



se2 



A collection of K different vectors {ci, C2, . . . , ck} in Z^ 
specifies a K dimensional subspace of 7i„ that is spanned 
by K graph-state basis {|Fc.)}f£i. 

For an example the graph state corresponding to the 
star graph S3 on 3 vertices with all edges weighted 1 as 
shown in Fig|T]^a) represents the GHZ state 



\9p^,)l^\j)2®\0p-j)3. (10) 



An edge with weight 1 will be represented by a black line 
and an edge with weight p — 1 will be represented by a red 
thick line as in Fig|2] In addition we will also indicate 
a vector in Z^ via colored vertices with white, black, 
blue, and red vertices representing weights 0,l,2,p— 1 
respectively. For example in the graph shown in Figjl|b) 
a vector (1, 0, 2) G Z^ is indicated via the colored vertices 
and therefore we have a graph-state base 
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III. CODING CLIQUES AND QECCS 

From the standard theory of the QECC we know that 
if a set of Pauli errors can be corrected then all the errors 
can also be corrected. The situation in non-binary case 
is the same, because the nonbinary error basis defined by 
& = {A'^Z*|s,t e Z^} is a nice basis [12j[16], i.e., the set 
E forms a basis for the operators acting on the n-qupit 
Hilbert space l-in. 

For a given weighted graph G — (V,r), from the defi- 
nition the vertex stabilizers of a graph state, we have 



t-s-r 



(12) 



That is to say every non-binary Pauli error acting on 
the graph state |r) can be equivalently replaced by some 
qupit phase flip errors, up to some phase factors. For 
convenience we refer that the vector t — s • F is covered 
by the error X^Z^ . Given an integer d we introduce a 
d-uncoverable set as 



- |t - s • F < |sup(s) Usup(t)| < dj , (: 



13) 



where we have denoted by sup(s) — {a £ V\sa ^ 0} the 
support of a vector s € and by |C| the number of the 
elements in C CV. That is to say is the set of all the 
d-uncoverable vectors, i.e., vectors cannot be covered by 
Pauli errors acting nontrivially on less then d qupits, in 
7^ . In addition we define the d-purity set as 



Sd = |s e Zp |sup(s) Usup(s • F)| < d| 



(14) 



It is obvious that if s e then the graph stabilizer Qs 
has a support less than d, i.e., act nontrivially on less 
than d qupits. 

A coding clique of a given graph G — {V,T) is a 
collection of K different vectors in Z^ that satisfy: 



i) e C^; 



K. 



ii) s ■ c = for all s e §d and every c e 

iii) c - c' e Dd for all c, c' e Cf . 

If the coding clique forms a group with respect to the 
addition modulo p, then it will be referred to as a coding 
group. In words, a coding clique is a collection of 
K d-uncoverable vectors in Z^ that is orthogonal to the 
vectors in Sd and the difference between any two vectors 
in is also d-uncoverable. We can build a super graph 
G with vertices being the vectors in Dd H Sj^, i.e., d- 
uncoverable vectors that are orthogonal to all the vectors 
in Sd, and two vertices in the super graph are connected 
by an edge iff their difference is also d-uncoverable. Then 
the coding clique Cf- is exactly a X-clique, a subset of 
the vertices which are pairwise connected, of the super 
graph G. This fact justifies the nomenclature clique. 

Theorem 1 Given a graph G = (V,F) and one of 
its coding clique C^, the subspace {G,K,d)p spanned 



by the graph-state basis {|rc)|c G C^} is an ((n, K, d))p 
code, which is a stabilizer code if is a group and a 
nonadditive code if is neither a group nor a subset of 
a coding group of G with d' > d. 

Proof. To prove that the subspace spanned by the 
graph-state basis {|Fc)|c £ Cf^} is an {{n,K,d))p code 
we have only to show that the condition |3D] 



(F,|£d|re') - f{£d)Sc 



(15) 



is fulfilled for all c, c' e and all the error £d = X^Z* 
that acts nontrivially on a number of qupits that is less 
than d, i.e., | sup(s) Usup(t)| < d. 

If the error is proportional to a stabilizer of the graph 
state, i.e, £d = f{s)Gs for some s £ Z^ and phase /(s), 
then s G §d- Because of condition ii) of the coding clique 
the error acts like a constant operator on the subspace 
(G, K, d) so that the condition Eq.( 15 ) is fulfilled with a 
f{£d) = /(■s) that is independent of c. If the error £d is 
neither one of the stabilizers of the graph state |F) nor the 
identity operator, i.e., 5d = X^Z* cx QsZ^~'''^ with t ^ 
s-T, then (Fcl^dlPc') oc (FlZ'^lF) with e = c'-c+t-s-T. 
By virtue of condition iii of the coding clique it is ensured 
that e ^ for all c, c' £ which gives rise to Eq.( 15 1 



with f{£d) = 0. Thus we have proved that {G,K,d)p is 
an ((n, K, d))p code. 

By definition, a nonbinary stabilizer code is the joint 
+ 1 eigenspace of an Abelian subgroup of generalized 
Pauli group P„ given in Eq.(j8]). The subspace (G, K, d)p 
spanned by the graph-state basis {|Fc)|c £ Cf^} is stabi- 
lized by an element S of P„ if and only if S is one of the 
stabilizers of the graph state Qg with s e § where 



{seZ^|s-c = 0,VceCf} 



(16) 



If 



is a coding group then it is obviously an 
Abelian group and can be generated by a set of vectors 
(ci,C2,...,Cfc) in Zjf with degrees [^i, /^2, • ■ • , Mfcli i-e-, 
the minimal positive number such that /^iCj — (mod p) 
{i = 1,2, ■ ■ ■ k). It is easy to see that all /i^'s divide p so 
that K — Y['i=i l^i divides p". Furthermore {ui = e*^'^/^') 



= m (17) 



cGt 



because £ Cd for the first equality and 



A; Mi-l 



E-^^-nE 



cGt 



K if s e §, 

if s ^ S, 



(18) 



(since fiiCi = (mod p)) for the last equality. That 
is to say if is a group then we can find exactly a 
number p"/K of stabilizers {Gsls £ S} whose joint -1-1 
eigenspace is exactly {G,K,d)p. Thus the code {G,K,d)p 
is a stabilizer code. 



If 



is not a group then we denote by C the group 
generated by C^, i.e., the smallest group that contains 
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Cf . Obviously the subspace Q spanned by {|rc)|c € C} 
is the joint +1 eigenspace of stabilizer {Gs\s € S} and 
|§| < p^/K. If the subspace Q can detect d — 1 or more 
errors then C is a coding group of G. If the subspace Q 
cannot detect d— 1 errors then any subset of the stabilizer 
cannot either. On the other hand every stabilizer code 
that contains {G,K,d)p must have a stabilizer that is a 
subset of the stabilizer of Q. Therefore if is not a 
subset of some coding group of G with d' > d then 
the code {G,K,d) is a nonadditive code. Q.E.D. 

Similar conclusions about the stabilizer codes appeared 
also in [38]. If all the generators of a coding group have 
the maximal degree p then corresponding stabilizer code 
can be denoted as [[n, /c, djjp. However there are cases 
where the generators of the coding group are not all of 
maximal degree. Then we have still a stabilizer code 
but the dimension of the code subspace may not be a 
power of p and we shall denote such a stabilizer code as 
{{n, ijLi-fi2-- ■ ■■tJ'k,d))p. A family of such kind of stabilizer 
codes will be provided in Sec. VI. 



IV. GRAPHICAL NONBINARY QECCS VIA 
NUMERICAL SEARCH 

According to Theorem 1 we can use the same system- 
atic algorithm developed for binary case in [2l] to do a 
systematic search for the non-binary quantum codes, i.e. 

i) To input a Zp- weighted graph G ~ {V,T) on n ver- 
tices; 

ii) To choose a distance d and compute the d-purity 
set Sd and the d-uncoverable set so that a super 
G can be built; 

iii) To find all the if -clique of the super graph G; 

iv) To output a {G,K,d)p, i.e., an {{n,K,d))p code 
that is spanned by the basis {|rc)|c G Cf-}. 

In practice, we have used the clique finding program 
diquer |39j to search for the chques for the super graph. 
Within our present computation power systematic search 
for graphical codes can be done up to p = 6 and n = 6 
and some tentative searches have been done for n — 8. 
In what follows a quantum code will be specified by a 
weighted graph together with a coding clique or the gen- 
erators of a coding group. Though the search for cliques 
are hard, the verifications of them are relative easy. 

A. The code [[3, 1,2]]3 

The first example is the stabilizer code [[3, 1, 2]]3, which 
is known and has been constructed, e.g., in |15j . We con- 
sider the star graph S3 on 3 vertices with two edges all 
weighted 1 as shown in Fig|lja). A coding group gener- 
ated by (1,0,2) e Z3 as shown in Figlllb) provides the 



code [[3, 1, 2]] 3, i.e., a 3-dimensional subspace spanned by 
the graph-state basis 

{Zi^Zf |53)|aeZ3}. (19) 

The stabilizer of this code is generated by {Q2, G1G3) with 
Qa being defined in ([6]) . It is easy to see that the stabilizer 
is equivalent to the stabilizer (<%'i23, Z123) appeared in [5] 
under a local unitary transformation. It is not difficult 
to see that we can have a [[3, 1, 2]]p for any odd p with 
the same graph and the same coding group generated by 
(1,0,-1) G Z3 with stabilizer (^a.^i^s)- 

B. The code ((3,3, 2))4 

For an odd number of qupits with p being even there is 
no code of distance 2 that saturates the Singleton bound 
so far. We have found a suboptimal code ((3,3,2))4 in- 
stead. For the equal weighted star graph on 3 vertices as 
shown in Fig|lja) we have found a coding clique: 

{(0,0,0), (1,0,2), (2,0,1)}, 

with corresponding graph-state basis shown in FigQa-c), 
which span the code subspace of a nonadditive ((3, 3, 2))4 
code. In fact, as we see later, we can construct a code 
{{3,p— 1, 2))p with even p. 

C. The code [[4, 2,2]]6 

The next example we concerned is also a 1-error detect- 
ing code [[4, 2, 2]] 6 which can be constructed from the star 
graph 5*4 on 4 vertices as shown in Figjljd). From this 
graph a 2 dimensional coding group can be found to be 
generated by vectors (1,-1,0,0) and (1,0,-1,0) in Zg. 
That is to say the code [[4,2,2]]6 is the 36-dimensional 
subspace spanned by the basis 

^Z-'^-'^Z'^Z^lSi) I a, 6 e Ze} . (20) 

whose stabilizer is generated by {GiG2G3,G4}- On the 
same graph we also find another coding clique as 

{ia + b,-a,^b,SaiSbi) eZl\a,beZe} (21) 

which can be stabilized by (^152^3) only thus we have 
a nonadditive ((4,36,2))g code. The nonadditive codes 
meeting the Singlet Bound is a very common situations 
in the nonbinary graphical codes, almost every stabilizer 
codes will have a clique set which are not a group, which 
means we can always construct a nonadditive code from 
a graphical nonbinary stabilizer code. 

D. The code [[5,1,3]]3 

The first example of 1-error-correcting code is the well- 
known [[5, 1,3]] 3 code, which is the only error-correcting 



FIG. 2: Twelve graph-state bases for the code ((6, 12,3))4. 
The thick red edge has a weight p — 1. 



codes for arbitrary dimension so far and was constructed 
in [17] and formulated with graph state in [30]. For 
the loop graph L5 on 5 vertices with all edges weighted 
1 we find a 1-dimensional coding group generated by 
(1,1,1,1,1) € Z| as indicated in FiglHe). It is clear 
that this coding group can be generalized to arbitrary 
dimension and we obtain for arbitrary p the [[5,l,3]]p 
code that is spanned by graph-state basis 



(22) 




(b) [[8,4,3]]3 



(c) [[8,2,4]]3 



FIG. 3: (a) Three generators of the coding group for the 
code [[7, 3,3]]3; (b) Four generators of the coding group for 
the code [[8,4, 3]]3; (c) All members of the coding groups for 
the code [[8,2,4]]3. 



E. The code ((5,4,3))4 

For dimension 4 we also found a nonadditive code 
((5,4, 3))4 that saturates the quantum Singleton bound. 
The graph we considered is still the loop graph L5 on 
5 vertices with all edges weighted 1. The coding clique 
contains the following 4 vectors in 

{ (0, 0, 0, 0, 0), (1, 1,1,1,1), (2, 3, 3, 2, 3), (3, 2, 2, 3, 2)} 

among which one basis is shown in FigQe) and two bases 
are shown in Figjll^f). Since there is a stabilizer of the 
code containing 4^ < 4^ elements only and it satuates 
the Singleton bound, the code is not a subcode of any 
1-error correcting code and is therefore nonadditive. In 
fact we can construct a nonadditive ((5,p,3))p for any 
p > 3 as will shown later. 



F. The code ((6, 12,3))4 

We also found a nonadditive code ((6, 12, 3))4 on the 
loop graph ig on 6 vertices with one edge weighted 3 
represented by a thick red edge as in Fig|2] The coding 
clique of Lg with 12 vectors reads 

(0 0) (0 1 1 2 1) (0 2 3 1 3) 

(1 1 3 1 0) (1 2 1 3 3 1) (1 3 3 3 2 3) 

(2 3 1 3) (2 1 3 2 1 1) (2 2 2 1 1 2) 

(2 3 2 3 0) (3 1 2 1) (3 1 2 2 2 2) 



which is represented in Figj2]with white, black, blue, and 
red vertices having weights 0,1,2,3 respectively. There are 
only 4 stabilizers can be found for this code, namely 

{I, G3GI, {QiQ2QAdbfQ?.Q&^ {QiQ2QAQbf Q^iGe,} ^ 

whose joint -1-1 eigenspace cannot be any 1-error correct- 
ing code because of the Singleton bound. Therefore the 
code is nonadditive. 



G. The codes [[7,3,3]]3 and [[8,4,3]]3 

We found also two families of optimal qutrit codes, 
namely the codes [[7, 3, 3]] 3 and [[8,4, 3]] 3, which can be 
constructed from the loop graphs and respectively. 
The generators of the coding groups are indicated in 
Figjsja) and Fig[3|b). For d = 3 the blue and red vertices 
coincide. However those generators of the coding groups 
are also vahd for the codes [[7, 3, 3]]p and [[8, 4, 3]]p for all 
odd p > 3 which will be discussed in details in the next 
section. 



H. The code [[8,2,4]]3 

The last example is the code [[8, 2, 4]] 3 which is found 
only recently |38j from the hype cube graph. The graph 
considered here is the equal-weighted wheel graph as 
shown in Figj3|c) and the corresponding graph state is 
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TABLE I: The codes for p = 3. 





d=2 


d=3 d=4 


n=o 


FFQ 1 oil 




n=4 


[[4,2,2]]3 




n=5 


[[5,3,21]3 


[[5,1,3]]3 


n=6 


[[6,4, 2]]3 


[[6,2,3]]3 


n=7 




[[7,3,3]]3 


n=8 




[[8,4,3]]3 [[8,2,4]]3 



denoted as jW^s)- The coding group is generated by 

(0,1,1,0,2,2,0,0) and (1,0,2,2,0,1,1,1) (23) 

and the graphical code (W8,2,4)3 is spanned by the 
graph-state basis 

{ZfZa'Zg^+^^Zf Z^'^+^'Z^^Zg^l W^s) |a, b e Z3} (24) 

The stabihzer of the code has the following set of gener- 
ators 

{G1G4, GaGe, G4G7, G4GS, GiG^Gz)- (25) 

To conclude this section we list all the stabilizer codes 
that saturate the Singleton bound, i.e., k + 2d < n + 2, 
in Table |I] It should be emphasized that the absence of 
some codes in Table [l] does not imply that these codes 
do not exist. To find some of these codes are beyond 
our computers' capacity and some of them will be con- 
structed in the following sections. Furthermore, we can 
find the codes with a length n > 6 via a random search- 
ing method, thus some of absent codes may still be found 
by using our systematic algorithm developed above. 



V. GRAPHICAL NONBINARY QECCS VIA 
ANALYTICAL CONSTRUCTIONS 



Because the clique finding problem is intrinsically an 
NP-complete problem, it is not plausible to rely on the 
numerical search for codes with larger length or higher 
dimension. Here we shall provide some analytical con- 
structions of good codes for higher dimension, which is 
based on the graphs and their coding cliques in lower di- 
mensions found via computer research. In practice, we 
start from a graphical code found for qutrit and then 
generalize the subspace to arbitrary dimension by adopt- 
ing the same graph and similar coding clique, and finally 
we prove that the subspace provides a code in arbitrary 
dimension. 



as shown in FigQa). The p — 1 dimensional subspace 
spanned by the basis 



(26) 



with < ^ < q — 1 and 0<j<q — 2isa. nonadditive 
code ((3,p - l,2))p. 

Proof — It is enough to demonstrate that a subset Cj^^ 
composed of 2g — 2 vectors of forms {1,0,21) and {q + 
j, 0, 2j + 1) satisfies all three conditions of coding clique. 
Since the 2-purity set is empty and I can assume value 
0, we have only to show that the condition iii of coding 
clique is satisfied. Any single qupit error can only cover 
those vectors of forms {a,b,a), (b,c,0), and (0,c, 6) with 
a,b,c £ lip being arbitrary. Considering the range of the 
I and j it is straightforward to show that all vectors in 
Cj""'^ and pairwise difference cannot be covered by single 
qubit error. Thus we have a code ((3,_p — l,2))p. The 
stabilizer of the code turns out to be generated hy Gb- 
Therefore it is a nonadditive code. 

Via a similar construction by Rains [23] we are able 
to construct the code {{2n + 3,p^"(p — 1), 2))p for even 
p. We consider the graph on 2n + 3 vertices composed 
of the star graph S3 as shown in FigQa) and the graph 
as shown in FigHib) and denote the corresponding 
graph state as j^s) ®\B2n)- Let us denote by {\v)} the 
basis in Eq.(26l for the ((3,p— 1,2))^ code constructed 
above then the code subspace is spanned by the basis 

(Z^Zc)^-i"^'(ZB)^"=i^--i|t;) ® Z^|B2n) (27) 

with s £ Zp" being arbitrary. We notice that the phase 
flips acting on |w) is a single qupit error on qupit B. 



B. The code [[2n, 2n - 2, 2]] 

We consider the graph on 2n vertices with all edges 
weighted 1 as shown in FigMb) and denote the corre- 
sponding graph state as \B2n)- The subspace spanned 
by the graph-state basis 



n-1 



n-1 



I[^2hl]lKi;U\B2n) (28) 



with aj, bj S Zp for j — 1,2, . . . n—1 is the code [[2n, 2n- 
2,2]] whose stabilizer is generated by 



. . . X2n-\Z2n, 
Z\X2Z3X^ . . . Z2n-\X2n- 



(29) 



A. The nonadditive code ((3,p— l,2))p with even p 

Description of the code — Suppose that p = 2q. We 
consider the star graph S3 labeled with V = {A, B, C} 



It is straightforward to see from the stabilizer that every 
single qupit error can be detected, i.e., not commute with 
at least one of two generators defined above. In compar- 
ison in [3S] the star graph has been used to construct the 
code [[n, n — 2, 2]]p. 
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FIG. 4: (a) The star graph 5*3 on 3 vertices; (b) The graph 
B2n', (c) The loop graph on 6 vertices with 1 edge weighted 
-1; (d) The loop graph Lt, (e) The loop graph Lg; (f) The 
wheel graph Ws on 8 vertices. 



C. The code [[2n + 3, 2n + 1, 2]]p with odd p 

Consider the graph on 2n + 3 composed of a subgraph 
S3 as shown in Figj4][a) and a subgraph B2n as shown in 
Figjij^b). On the first 3 qupits a [[3, 1, 2]] code can be con- 
structed for odd p with the stabihzer given by {QbiGaQc) 
as shown previously. Then we obtain the stabihzer of the 
code [[2n + 3, 2n + 1, 2]]p as 



ZaXbZcXiZ2X3Z4^ . . 



(30) 



TABLE II: The stabihzer of the code [[6, 2, 3]]p 





1 


2 


3 


4 


5 


6 




X 


z 


z 


X 


Z 


z 


G3G6 


Z 


z 


X 


z 


Z 


X 


Q1Q2QZ 


zx 


xz'^ 


zx 


z 


X 


z 




X 


z 


zx 


xz'' 


zx 


z 



not a multiple of 3 and 3g^ if p is not 3 but can be divided 
by 3. In either case the number of stabilizers is less than 
p"* if p > 3 and since the code saturates the Singleton 



bound we conclude that the code Eq.(32 1 is nonadditive 



E. The code [[6,2,3]]p 

Description of the code — Consider the weighted loop 
graph on 6 vertices as shown in FigQc) and the corre- 
sponding graph state \Le). The p^-dimensional subspace 
spanned by the graph-state basis 



,beZp} (33) 



is a [[6, 2, 3]]p code for all odd p > 2. A set of the gener- 
ators of its stabilizer is listed in Table |lTl 

Proof — Because of Theorem 1 we have only to prove 
that the subset of vectors 



D. The nonadditive code ((5,p, 3))p with p > 3 

We consider the loop graph L5 with all edges weighted 
1 and the subset C| C Zp composed of p vectors including 
{jijiiijii) with j ^ 2,-1 and two additional vectors 
a = (2,-1,-1,2,-1) and b = (-1,2,2,-1,2). All of 
these vectors have full support and obviously cannot be 
covered by any single qupit error which is able to cover 
at most 3 vertices. All vectors (ci,C2,--- ,05) with full 
support that can be covered by a 2-qupit error has the 
property 



(31) 



It is obvious that all vectors in C3 do not have such a 
property. Furthermore the difference of any pair of vec- 
tors in C3 has full support and do not have the property 
Eq.(31 1 either. Therefore C3 is a coding clique since con- 
ditions i and ii are trivially satisfied. To find all the 
stabilizer of the code constructed from the coding clique 
i.e., subspace spanned by the graph-state basis 



-3' 



{(Z1Z2 232425)^1X5), 2" lis), 2^|L5)}, 



(32) 



we have to solve equations in Eq.(16). I t tu rns out that 
the total number of the solutions to Eq.(16) is q'^ if p is 



{c = (a, a + b, b, —a, a ~ b, b)\a, b £ Zp} (34) 

satisfies all three conditions of the coding clique. Con- 
ditions i) and ii) are obviously satisfied since a, b can 
assume value zero and the 3-purity set is empty for the 
loop graph for any dimension p. To prove condition iii) 
we have to show that c cannot be covered by single and 
2-qupit errors for all a,b G Zp. Obviously c cannot be 
covered by any single vertex error, we consider only 2- 
qupit errors in what follows. 

i) = 0, a 7^ with c = (a, a, 0, —a, a, 0). If a 2-qupit 
error is supported on vertices {1, 2, 4, 5} then it can 
only takes form 



x^z^x.-zl 



■yn -yji '-^t 



-^^3^6, 



which covers (s, n, 0, n, t, 0), (n, s, 0, t, — n, 0) and 
(— s, t, 0, i, s, 0) respectively with n,s,t G Zp being 
arbitrary. All these 3 types of vectors cannot be 
identified with c provided p > 2 and odd. For ex- 
ample if (s, n, 0, n, i, 0) = c then 2o = which is 
impossible for odd p. 

ii) a = 0, 6 7^ with c = (0, 6, 5, 0, —6, b). In this case 
it can be proved in exactly the same manner as the 
first case that c cannot be covered by 2-qupit error. 
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iii) a = b ^ with c = (a, 2a, a, —a, 0, a). Any 2-qupit 
error that covers {1,2,3,4,6} can only take form 
X^Z^X^Zt^ or X^Z^X^, which covers 

(s, m + n, n, 0, — m) or (to, s, to, n, 0, n) 

respectively with to, n, s,t € Zp being arbitrary. To 
cover c we have to require 2a = —2a or a = —a 
which are impossible for odd p. 

iv) a = -b ^ with c = (-5, 0, 6, 6, -26, &). Any 2- 
qupit error that covers 5 vertices {1,3,4,5,6} can 
only be X:^'X^'ZI or XY'ZiXg'Z*^, which covers 

(to, 0, TO, n, s, n) or (— n, 0, to, s, to + n, t) 

respectively with to, n, s,t € Zp being arbitrary. To 
cover c we have to require b = —b or 26 = —26 
which is impossible for odd p. 

v) a ^ ±6, a, 6 ^ with c = (a, a + 6, 6, —a, a — 6, 6). 
To cover all the 6 vertices a 2-qupit error can only 
takes form 

which covers {s^m,n,t,n, —m), {m, s,m,n,t,n), 
and (—n,m,s,m,n,t) respectively. However all 
these vectors can not be identified with c if p is 
odd. For example if (s, m, n, t, n, —to) = c then 
2a = which is impossible for odd p. 

In summary, for any a,b ^ "Lp the vector c cannot be 
covered by any single or 2-qupit errors so that is a coding 
clique and corresponding subspace is the [[6, 2, 3]]p code 
for all odd p> 2. 



F. The code [[7, 3, 3]]p 

Description of the code — Consider the equal weighted 
loop graph on 7 vertices as shown in Figj4][d) and the 
corresponding graph state \L-j). The p^-dimensional sub- 
space spanned by the graph-state basis 

{Zl+''+-Z^ZlZlzr-Z^-Z',\Lj)\a, b,ceZp} (35) 

is a [[7, 3, 3]]p code for all odd p > 2. A set of the gener- 
ators of its stabilizer is listed in Table IIIII 



TABLE III: The stabilizer of the code [[7,3,3]]p 





1 


2 


3 


4 


5 


6 


7 




I 


Z 


X 


Z 


Z 


X 


Z 


G4G7 


Z 


X 


z 


X 


Z 


z 


X 


G2G3G5 


z 


XZ 


zx 


Z' 


X 


z 


X 


G1G2G3G4 


xz 


X 


xz"^ 


ZX 


z 


X 


z 



Proof — Because of Theorem 1, we have to show that 
the subset of Zp defined as 

{c— (a + 6 + c, a, c, 6, a — c, — c, 6) |a, 6, c £ Zp} (36) 

is a coding clique. Since the 3-purity set is empty for the 
loop graphs with any dimension p, we need only to show 
that any nonzero c cannot be covered by any single or 
2-qupit error. 

It is not difficult to see that c cannot be covered by any 
single qupit error, which covers a vectors with 4 consec- 
utive components being zero. Because of the symmetry 
of the loop graph, all 2-qupit errors can be classified into 
3 types Ti2,Ti3, and T14 where Tij denotes two errors 
occur on qupits n + I and n + j with n G Z7. Each type 
of errors covers the vector (ci, C2, . . . , cy) with properties 



T12 

Ti4 



Cn+l = Cn+2 — 0; 

+ c„+3 = c„+i and c„_2 = c„_3 = 0; 
= c„+i, c„_2 = 0, and c„_3 = c„+2- 



It can be easily checked that as given in Eq.(36 1 for arbi- 
trary a, 6, c e Zp does not belong to all these 3 types of 
vectors when p is odd. For example a T12 type of error 
covers vector with 3 consecutive components being zero, 
which is impossible to be identified with a nonzero c. 
That is to say c cannot be covered by any 2-qupit errors 
and therefore we have a coding chque and a [[7,3,3]]p 
code for all odd p. 



G. The code [[8,4,3]]p 

Description of the code — Consider the equal weighted 
loop graph on 8 vertices as shown in Figj4](e) and the 
corresponding graph state \Ls). The p^-dimensional sub- 
space spanned by the graph-state bases 



^1^2 



c ^e— c 
^3 



^4 ^5 ^6 



ZfZg^-^+^lLs) (37) 



with a, b,c,e E Zp is a [[8, 4, 3]]p code for all odd p > 2. A 
set of the generators of its stabilizer is listed in Table IV 



TABLE IV: The stabilizer of the code [[8, 4, 3]]p 





1 


2 


3 


4 


5 


6 


7 


8 


GiGaGiGa 


Xz 


X 


XZ 


ZX 


Z 


X 


Z 


ZX 


GiGiGsGc) 


z'' 


X'^Z 


Z^X 


Z' 


Xz 


ZX 


Z 


X 


/-» 2 /-» 2 /-» 2 /-» 

t/lt/4t/5t/7 


x^ 


z' 


z' 


X'^Z'^ 


Z'^X^ 


z 


X 


z' 


G1G2G3G4G5 


ZX^ 


ZX 


X 


xz'^ 


XZ 


z 


X 


z^ 



Proof — Because of Theorem 1, we have to show that 
the subset of Zp composed of vectors of form 

c = (e, 6— c, e— c, e— a, a-|-6, a—b+c+e, 26, a—c+e) (38) 
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which with a, b,c,e G Zp is a coding clique of the loop 
graph. Because of the symmetry of the loop graph, all 
2-qupit errors can be classified into 4 types Ti2,Ti3, 
Ti4, and T15. Each type of errors covers the vector 
(ci, C2, . . . , cg) with properties (n S Zg) 



T12 

Ti3 
Ti4 



Cn+l = Cn+2 = Cn+3 = Cn+i = 0; 

Cn+2 = Cn+3 = c,i+4 = and C„+i + C„_3 = c„_i; 
Cn+2 = Cn+3 = 0, c„+i = c„_i, and c„_2 = c„_4; 

Cn = C„+4 = 0, Cn+l = C„+3, and Cn-1 = C„_3. 



As an example we consider a T12 type of errors which cov- 
ers a vector with 4 consecutive components being zero. 
A nonzero vector c as defined in Eq. ( 38 1 is not possible 
to have such a property, e.g., the last 4 components be- 
ing zero. In the same manner it can be checked that c 
does not belong to all those 4 types of vectors above for 
arbitrary a, b,c,e € Zp when p is odd. Since all single 
qupit errors cover vectors with the same property as in 
type T12, we conclude that c cannot be covered by any 
single or 2-qupit errors so that we have a coding clique 
and therefore a [[8,4, 3]]p code for all odd p. 



H. The code 



.2 AT 



Description of the code — Consider the equal weighted 
wheel graph on 8 vertices as shown in Fig. 3(d) and the 
corresponding graph state \ Ws)- The p'^-dimensional sub- 
space spanned by the graph-state basis 

{z1Z^Zr''2fzl^Z',-^ZllZ',\T)\a, b e Zp} (39) 

is a [[8,2,4]]p code for all odd p > 2. A set of the gener- 
ators of its stabilizer is listed in Table Ivl 



TABLE V: The stabilizer of the code [[8, 2, 4]]p 





1 


2 


3 


4 


5 


6 


7 


8 


G1G7 


X 


z 


z 


X 


Z 


z 


X 


X 


QiQs 


XZ 


z 


X 


z 


Z 


X 


z 


ZX 


G1G4 


X-' 


z' 


z 


X 


Z 


X 


X 


z 


Q2G5 


z 


x^ 


z' 


z 


X 


z 


X 


X 




X 


z^ 


X 


z 


z 


X 


z' 


X 


G1G2G3 


zx 


xz'^ 


zx 


z 


z 


z 


z 


z 



Proof — Because of Theorem 1 we have to show that 
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the subset of Zp defined as 



{c ^ (6, a,a-b, 2b, 2a, b~a,b, b) \a, belp}. (40) 

is a coding clique of the wheel graph. We have only to 
prove condition iii of the coding clique since the 4-purity 
set is empty. 



A single qupit error on vertex n can cover a vector with 
property c„±2 = c„±3 = which is impossible for c de- 
fined above. Because of symmetry all 2-qupit errors can 
be classified into 4 types T12, T13, T14, and T15 with each 
type of errors covering the vector (ci, C2, . . . , cg) with 
properties (n G Zg) 



T12: 


Cn- 


f 1 = Cn- 


f4, C„+2 


= c„+5 = 0, and c„+3 = c„+6 


T13: 


Cn 


= Cn+3 


and Cn+ 


1 = c„+4; 


T14: 


Cn 


= Cn+3 


and Cn+ 


1 = c„+2 = 0; 


T15: 


C-n 


C-n+Q 


Cn+2 — 


c„+4 and Cn+l = c„+5 = 0. 



It can be checked in a straightforward manner that none 
of the above equalities can be satisfied by c ad defined in 
Eq.(40). For example we consider error of type T12 with 
n = 0. In this case we have constraints b — 2b, a = 2a, 
and a — b = b — a which are impossible for nonzero c when 
p is odd. This is true for all n G Zg. Thus any 2-qupit 
error cannot cover c. 

All 3-qupit errors can be classified into 7 types 
iri23,Tri24, Tri25, T126, T127, T135, and Ti36 with each 
type of errors covering vectors with the following prop- 
erties (71 G Zg) 



'^123: 


Cn- 


f 2 = Cn 


^5 and 


Cn-2 


— Cn 


-5; 


1424: 


Cn 


= Cn+3 


and Cn 


+2 = 


0; 




T125: 


Cn 


+ Cn+l 


= C„+3 


and 


Cn+4 


= 0; 


'ri26: 


Cn 


+ Cn-1 


= C„_3 


and 


C„-4 


= 0; 


'iri27: 


Cn 


= Cn-3 


and Cn 


-2 — 


0; 




T135: 


Cn 


+ c„+i 


= C„_|_3 


and 


Cn + 


C„-l 


Ti36: 


Cn 


— Cn+2 


+ C„+3 


— Cn 


-2 + 


Cn-3 



^n — 35 



It can be checked in a tedious but straightforward manner 
that for every n G Zg none of those equalities above can 



be satisfied by c as defined in Eq.(40). For an example 
we consider the 3-qupit error of type T126 with n = 1. 
In this case we have b + b = b — a and 2a = which are 
impossible for nonzero c. 



To summarize, the vector c defined in Eq.(40l for all 
a, 6 G Zp is 4-uncoverable. Thus we have proved that 
the subspace defined in is a [[8,4, 2]]p code for all odd 
dimension p. 



VI. CODES FROM COMPOSITE SYSTEMS 

Consider a system with pq levels whose computational 
bases are denoted by {[Opglflo ^- We can also regard this 
system as a composite system of a p-level system and a q- 
level system, whose computational bases are denoted as 
{|s)p}s=o ^^'^ {\i)q}t=o respectively. If there are n copies 
of pq-level systems, we also have n copies of p-level and 
q-level systems. On the other hand if we have two groups 
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of p-level and g-level systems we can also obtain n copies 
of pg-level system by pairing up one p-level system and 
one g-level system to made up a composite system. 

Given a Zp-weighted graph {V, Tp) and a Z^-weighted 
graph (y, Tq) on the same vertex set V and corresponding 

coding cliques and C^, the subspace spanned by the 
basis 



c e 



c e 



'4} (41) 



is an {{n, KK,d))pq code with n = \V\. This is because 
all the direct products of Pauli errors of the p-level sys- 
tems and g-level systems form a nice error basis for the 
pq-level system. That is to say, via a direct product of 
two graphical codes (G, K, d)p and (G, K, d)q we can con- 
struct a code {{n, KK , d))pq, which is however not nec- 
essarily to be another graphical code. As will be shown 
below this construction will yield a graphical code of a 
higher dimension if p and q are coprime, in which case 
there exist two integers a and (3 such that 



ap + (3q— 1. 



(42) 



Given a Zp^-weighted graph (V, Tpq) we can build a Zp- 
weighted graph {V,Tp) and a Z^-weighted graph {V,Tq) 
whose adjacency matrices are given by 



pTpq = Tq (mod q), qVpq = Tp (mod p). 



(43) 



On the other hand, given a Zp-weighted graph (VjFp) 
and a Z^-weighted graph {V,Vq) on the same vertex set 
V, we can also build a Zp^-weighted graph iV^Tpq) with 
adjacency matrix given by 



r 



pq 



2t^ 

pa Iq 



ql3^Tp (mod pq). 



(44) 



By relabeling of the bases of a pq-level system according 
to \s)p ® \t)q i-^ \pt -\- <is)pq, we cau define an isometry 
between n pq-level systems and n pairs of p-level subsys- 
tems and g-level subsystems as 



n 



E E b*- 



qs) 



19' 



(45) 



which is only possible when p and q are coprime, we have 



\Tpq)^n\Tq). 



(46) 



Accordingly the bit flips and phase flips are related to 
each other via 



Xpq — Ti-iX^ 



x")n\ Zpq = niZp®Zq)nl (47) 



Theorem 2 If p,q are coprime then the graph state 
on Zpq-wcightcd graph is in a one-to-one correspondence 
with the direct product of two graph states on a Zp- 
weighted graph and a Zg- weighted graph whose adja- 
cency matrices are related via Eqs.( 43|44 ). 

According to the fundament theorem of arithmetics 
any integer p can be expressed as p = p"^P2 ^ . . . for 




(b) 



+ 
a. 



+ 




4-- 



FIG. 5: An example of the composite system codes, (a) The 
inner graph is a Z2-weighted graph and the outer graph is a 
Zp-weighted graph, (b) The Z2p-weighted graph Ha where 
5=(p+l)V2. 



some primes pi {i — 1, . . . , L). Therefore a graph state 
on the subsystems of dimension p will be in a one-to- 
one correspondence of a direct product of graph states 
on the subsystems of dimension p"' . As a consequence 
we have only to consider the graph states of systems of a 
dimension that is prime or a power of prime. So far the 
deflnition of graph state is unique in the case of prime 
dimension and there is a second definition of graph states 
[55] in the case of prime power dimension. 

As an example we consider the case where q = 2 and 
there 6 copies of 2p-level composite systems and accord- 
ingly 6 pairs of qubits and qupits. A direct product of a 
graphical code [[6, 1, 3]] 2 on qubits constructed from the 
loop graph with a leaf [57] as shown in Fig. 5(a) and the 
graphical code [[6,2,3]]p arising from the graph Lg as 
constructed in Sec.V(E) will result a code ((6,2 -^^,3)). 
It is in fact a stabilizer code and we should note that the 
dimension of the coding subspace is not a power of 2p. 

Via the relabeling of the bases of the 2p-level systems 
as specified TZ in Eq. ( 45 1 we obtain a Z2p- weighted graph 
Hq as shown in Fig. 5(B). From this graph we consider a 
subset 



of 2p^ vectors in Z2p of form 



(p -I- l)(a, a + b,b, —a, a — b,b) + p{c, c, c, c, c, 0) (48) 

with a^b Tjp and c e Z2. The subspace spanned by the 
graph-state basis {Z^\Hq)\c e C} is exactly the stabilizer 
code ((6, 2p^, 3)) constructed from composite systems. A 
set of the generators of its stabilizer is listed in Table VI. 

Because the code [[6,2,3]]2 does not exist we cannot 
remove any stabilizers from the lower half part of the 
above table to obtain a (2p) ^-dimensional code space, 
therefore the code ((6, 2p^, 3))2p is not a subcode of a 
[[6, 2, 3]]2p, if there is any. 



VII. DISCUSSION 

We have generalized the graphical construction JTT of 
QECCs to the nonbinary case to find both additive and 
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TABLE VI: The stabilizer of the code ((6, 2p^ 3))2p. 





1 


2 


3 


4 


5 


6 
















QlQl 












X^ 


/-^2 /-^2 










X 




QiQ^Qh 


X 








Z^+^X^ 




Ql 


X 


J 


X 


J 


z^ 




QlQl 




Z^X^ 


Z^ 


X 


z^ 


X 


QlQl 


X^ 


X 


x^ 


Z^ 


Z'' 


X 


QlQl 


X^ 


z^ 


z^ 


x^ 


X 


X 


QlQl 


X'PZf 


Zf 


X 


Zf 


ZPX" 


zp 



nonadditive codes based on nonbinary graph states. The 
advantages of our graphical approach lies in the fact that 
we are able to construct codes on physical systems of ar- 
bitrary dimension, prime or nonprime, to construct both 
additive or nonadditive codes, pure or impure. In ad- 
dition the basis for all codes are explicitly constructed. 
In principle for prime dimension our method exhausts 
all the stabilizer codes, which can be demonstrated in 
exactly the same manner as in binary case |27j . 

Via numerical search we have found many opti- 
mal codes including two optimal codes [[8,4, 3]]3 and 
[[8, 2, 4]]3 which have been found in [38 on a hyper cube 
graph. Since the clique finding problem for non-binary 
case is even more difficult and is intrinsically a NP- 
complete problem, our computational result is mainly 



limited in codes with n < 8 and small distance d < 4. 
With the help of the codes found numerically we also 
manage to construct analytically some families of opti- 
mal stabilizer codes that saturating the Singleton bound 
such as [[6,2,3]]p, [[7,3,3]]p, [[8,4,3]]p and [[8,2,4]]^ for 
any odd p. We have also explicitly constructed the code 
[[n, n — 2, 2]]p except the case of even p and odd n. There 
exist stabilizer codes whose code subspace is not a power 
of the dimension of physical systems such as the code 
((6, 2 -p^, 3))2p which is constructed via a composite sys- 
tem approach. Furthermore we have constructed a family 
of nonadditive codes ((3,p — l,2))p with even p, a non- 
additive optimal code ((5,p, 3))p for all p > 3. 

Finally, we briefly address some questions which are 
still open. Although our method can also be used to find 
additive and nonadditive codes for non-prime dimension, 
we can not exhaust all the non-prime stabilizer codes be- 
cause of our definition of graph states. Because of The- 
orem 2 we have to consider the graph states on system 
of prime and prime power dimension. There are at least 
two different definitions of the graph state on n subsys- 
tems of a dimension p = q'" being a power of prime: one 
is as defined in this paper and in [38^ and the other one 
is as defined in [33] . A different definition of graph states 
may result in some new families of codes. 
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